Where is the spectral weight in magnetic neutron scattering in the cuprates? 
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We present estimates in the Hubbard and Heisenberg models for the spectral weight in magnetic 
neutron scattering experiments on the cuprates. With the aid of spin-wave theory and the time 
dependent Gutzwiller approximation we discuss how the spectral weight is distributed among the 
different channels and between high and low energies. In addition to the well known total moment 
sum rule we discuss sum rules for each component of the dynamical structure factor tensor which 
are peculiar for spin 1/2 systems. The various factors that reduce the spectral weight at the rel- 
evant energies are singled out and analyzed like: shielding factors, weight at electronic energies, 
multimagnon process etc. Although about 10% ~ 15% of the naively expected weight is detected 
in experiments after consideration of these factors the missing weight is within the experimental 
uncertainties. A large fraction of the spectral weight is hard to detect with present experimental 
conditions. 
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I. INTRODUCTION 

Magnetic neutron scattering (MNS) in high temper- 
ature superconducting cuprates usually detects about 
10% ~ 15% of the spectral weight dictated by a naive 
application of sum rules. For example the total weight 
in a wide range of energy and momentum in a recent 
experimenti in La2-2:Ba2:Cu04, with x — 0.125, is ~ 
0.22/x^ whereas in the insulating phase the naive expec- 
tation from sum rules is that one should find 2/i^. It is 
usually argued that this value should be corrected for the 
hole destruction of moments by a 1 — x factor which still 
leaves a large fraction of spectral weight undetected. 

This rises various problems in the interpretation of 
MNS. For example it has been argued that the average 
of the dynamical susceptibility weighted by the Fourier 
transform of the magnetic interaction can be used to es- 
timate the energy involved in magnetic pairing and its 
temperature dependenceiSii^i^ Clearly to obtain an ab- 
solute estimate the spectral weight problem needs to be 
sorted out first. Furthermore modeling the dynamical 
structure factor probed by MNS becomes rather prob- 
lematic since sensible theoretical models do satisfy sum 
rules. Indeed any theoretical claim of intensity agree- 
ment with the measured dynamical structure factor in 
absolute units needs to explain how the sum rule is sat- 
isfied or why it is violated. This is even more stringent 
in spin only models for which neither 1 — x factors nor 
shielding corrections apply. 

The purpose of this work is to explain this apparent 
discrepancy. We provide theoretical estimates of the vari- 
ous factors which correct the sum rule and estimate what 
fraction of the spectral weight is accessible to present 



day experimental conditions. Theoretical estimates are 
performed in the antiferromagnetic (AFM) phase using 
the Heisenberg and the Hubbard model combining spin- 
wave theory, numerical results and the time dependent 
Gutzwiller approximation (TDGA)£i2iS and in the doped 
phase in the Hubbard model within the TDGA. Apart 
from the mentioned I — x factor we discuss the so called 
"shielding factors" due to an incomplete formation of 
magnetic moments. We estimate the spectral weight loss 
to electronic transitions at energies too high to be de- 
tectable by present day inelastic magnetic neutron scat- 
tering experiments and also the weight in multimagnon 
processes which is either at too high energies or is so 
broad in energy and momentum that it is not detectable 
in unpolarized neutron scattering experiments. After 
consideration of all these factors we arrive to the conclu- 
sion that within the experimental uncertainties the sum 
rule is not violated (which is reassuring) on the other 
hand a major fraction of the spectral weight is very hard 
to detect with present experimental conditions. 

The outline of the paper is at follows. In Sec. ITTl we 
shortly review the theory of magnetic neutron scatter- 
ing and the relevant sum rules to fix notations. This 
section has also a pedagogical character. Apart from 
the well known total moment sum rule we discuss sum 
rules for each component of the dynamical structure fac- 
tor tensor which, to the best of our knowledge, have not 
been applied in the present context. We also highlight 
some simple experimental facts that are usually assumed 
as granted in experimental works, like domain averages 
(Sec. Ill but often overlooked in theoretical works. In 
Sec, mil we discuss the spectral weight distribution in the 
undoped case and in Sec. II VI we discuss the doped case. 
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We conclude in Sec. |V] 



II. MAGNETIC NEUTRON SCATTERING 

We start with a short review of magnetic neutron scat- 
tering to fix notations and discuss the sum rules that are 
relevant to our problem. 

The magnetic neutron scattering cross section is given 



dfldE' k \2nB J 

X ^((5„;3-g„g;3)5"'3(g,,,) (1) 
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where e"^^'^''^ is the Debye- Waller factor, q = k—k' {q = 
q_/\q\), k {k') is the initial (final) wave vector of neutrons 
and {'^r^/2^BY = 72.65 x 10~'^barn//i^. The magnetic 
form factor is given by F{q) = J (ire*'' ''|(/)(r)|^ where 
(j){r) is a Wannier orbital and we defined the dynamical 
structure factor tensor, 

^"'(9'-)-^ (2) 
X e-P^^{^i\S1^\v){v\S^^\^Ji)5{hu; -E, + E^l 



where g is the Lande g- factor. For free electrons g = 
2.0023. In a solid a different value may be appropri- 
ate which may also depend on direction. For example 
Ref. 11 quotes g = 2.08 in the plane and g — 2.36 per- 
pendicular to the plane for a typical cuprate. For sim- 
plicity we take an isotropic g unless otherwise specified. 
S'q is the Fourier transform of the a component of the 
spin operator and Z is the partition function. The dy- 
namical structure factor S"''^{q, oj) obeys detailed balance 



A. Sum rules 

We will discuss sum rules for effective models of the 
magnetic dynamics. Because effective models restrict the 
Hilbert space sum rules turn out to be model dependent 
and therefore should be applied with care on modeling a 
real system, as discussed below. Our considerations are 
based on the two most popular models in this context, 
namely the Heisenberg and the Hubbard model, respec- 
tively. 

The so called total moment sum rule is usually formu- 
lated within a Heisenberg model with spin S and reads: 

= ^ E d{ncj)S"^{q, oj) = {gfiB fSiS + 1) 

(3) 

This applies to a system where magnetic ions have one or 
several partially filled orbitals (for example a rare earth 



ion with a partially filled /-shell) and well formed mag- 
netic moments i.e. when double occupancy of a given 
orbital is negligibly. 

Within this work we will restrict to systems where ions 
have only one partially filled orbital and the system can 
be model with a one-band Hubbard model: 



H 



t' 



E/ ^i<J^3<y 



Uy^m^n^i. (4) 



Here cj^. (qo-) destroys (creates) an electron with spin a 
at site i, and riia- = c\^Cia- U is the on-site Hubbard 
repulsion and both nearest (~ t) and next-nearest (~ 
t') neighbor hopping has been included. Most of our 
considerations apply also to other models where ions have 
only one partially filled orbital per atom, like the usual 
three-band Hubbard model for cuprates with Cu d and 
O p orbitals. 

We define for later use the spin autocorrelation func- 
tion and the zeroth moment of the diagonal components 
of the dynamical structure factor as: 



d(nuj)S"°'{L0) 



(5) 
(6) 



When one has one orbital per site (i.e. a spin 1/2 
system) more stringent sum rules than Eq. (PJ apply. 
Indeed each component of the dynamical structure fac- 
tor satisfies a separate sum rule. For example from 
Eqs. one finds: 



{n - 2D) 



(7) 



where S*^ = e-^i '-^Sf with Sf = (n,| -ni|)/2 and we 
used that nf^ = ma. We also defined the thermal ((...)) 
and spacial averages of the orbital occupancy: 



and double occupancy: 

D = ^E^'^^T'^^i)- 

i 

Analogous proofs hold for the other components: 

M^ = fil{n~2D) (8) 

where we took g = 2. (We will occasionally restore g 
below when convenient for clarity). Eq. also follows 
from the fact that the quantization axis in Eq. Q is 
arbitrary. 
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Eq. Q is valid in the Heisenberg model where D = 
and S is arbitrary. Eq. ||SJ) in contrast is valid for S* = 1/2 
systems but without restriction in D. 

The factor {n — 2D) in Eq. ||SJ) is the probability to find 
an atom singly occupied and reflects the fact that dou- 
bly occupied or empty atoms do not produce magnetic 
scattering of neutrons reducing the total scattering cross 
section. This is some times called "shielding factor" . 

One can also prove that the total weight of the off- 
diagonal components of 5"'^ adds to zero. First no- 
tice that only the symmetric part of S°'^ contributes to 
Eq. (|TJ: 

J2iSc.f3-q»qf3)S'^''{q,u;) 

a(3 

Using the Lehmann representation Eq. (0) one can show 
that: 

1 roo 

q •> 

= S^pfil{n-2D) (9) 

which gives us a sum rule for each component of 

We are not aware of references quoting the sum rule 
Eq. although the result is so simple that we doubt it 
is original. 

If S^gf = J2i is ^ good quantum number S°'^{q,Lj) 
becomes diagonal with S^^{q,uj) = Syy{q,u!)m^ In this 
case we can write the cross section as 

X [(1 - + (1 +g)S-(9,i^)] ^ 

We will restrict to systems where this expression applies. 



B. Ordered states and Bragg scattering 

In the presence of long-range magnetic order the sys- 
tem shows magnetic elastic scattering. We will consider 
phases in which spin rotational invariance is broken with 
order along the z axis (i.e. stripes, AFM, etc.). We as- 
sume a magnetic unit cell with atoms (A^*^ = 2 for 
the AFM state at half- filling) at positions R + 6i, with R 
the cell position and 6i the position of the atoms within 
the ceh {i = 1,...,N^'). The vectors R fo rm a Bravais 
lattice. For such a magnetic structure we have: 

{Sq) = N^Sq^Q^.TUQ^, 
Qm 



where the sum is over the magnetic reciprocal basis vec- 
tors and we defined: 

"^QM^^^Ee^"^"'''^- (11) 

" i 

and the local site-dependent magnetization rrii — (n^i — 
n^)/2. 

It is convenient to define the fluctuation operator: 

63^ ^ SI - 

With these definitions the longitudinal structure factor 
can be put as: 

Qm 

+ Y,M^W\''^i^^-E-+Eo)] (12) 

where for simplicity we set T = and S{q) is Kronecker's 
S whereas S{?iuj) is Dirac's S. The first term in the brack- 
ets describes Bragg peaks. The weight of the peaks is 
given by the square of the Fourier transform of the mag- 
netization inside the magnetic unit cell Eq. Hll|l . The 
second term describes inelastic scattering. The inelastic 
part (only!) of the dynamical structure factor is related 
to the dynamical susceptibility via the fluctuation dissi- 
pation theorem. 

C. Domain average 

For the case discussed above of a diagonal struc- 
ture factor tensor, the cross section involves the factors 
(1 — q^). These are polarization factors for scattering 
with unpolarized neutrons and are rooted in the dipolar 
interaction between neutrons and the electron spin. In 
experimental works often the polarization factors (1 — 9^) 
are included in the definition of the dynamical structure 
factor and an average over the orientation of domains, 
{...)dom, is done: 

S^fHq,cu) = ^((l-g2)),_^-(g,^) (13) 

a 

For a paramagnet all directions are equivalent, S°'°'{q, uj) 
does not depend of a and ((1 — q^))dom = 2/3. Ordered 
systems will be characterized by an order parameter that 
is a vector like the staggered magnetization. In general 
a real sample will consist of domains with different ori- 
entations of the order parameter. For a distribution of 
orientations that is completely isotropic in spin space 
((l-€))do™ = 2/3. 

Consider now scattering from a quasi two-dimensional 
(2d) system with the c direction defined perpendicular to 
the plane and the a and 6 directions in the plane. In quasi 
2d-systems the most common experimental configuration 
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is that the planes are perpendicular to the incident neu- 
tron beam and, depending on the energy, the component 
of q perpendicular to the plane, Qc may be larger than 
the components in the plane. In the extreme case that 
Qc >> qa,<jb one can put [c.f. Eq. (|13() ]: 



(14) 



This is valid regardless of the domain distribution. In 
order to evaluate this expression it is convenient to use a 
domain-dependent reference system in which the z axis 
follows the ordered moment of the domain. In the case 
of an isotropic distribution of domains one recovers (1 — 
ql))dora = 2/3. That is: 



S^ff{q,oj)^-[S-^{q,Oj) 



-Syy{q,Lo)] + -S^\q,Lo) (15) 



where we have grouped the transverse contribution. On 
the other hand, one usually deals with systems that have 
easy planes or easy axes. In cuprates for example the 
ordered moment is usually in the Cu-0 plane or close 
to that plane. For a distribution of the ordered moment 
which is isotropic within the plane and in the above scat- 
tering geometry 



-Syy{q,u;)] + S^%q,uj). (16) 



The 1/2 factor is due to the fact that one transverse mode 
is perpendicular to the plane and becomes silent in this 
configuration. 

For reasons that will be clear below real experiments 
detect mainly the transverse structure factor. Comparing 
Eqs. ifTsl) . ((TB)l we see that in order to accurately estimate 
spectral weights from experiment precise information is 
needed on the domain distribution since the transverse 
component appears with different weight. 

The condition qc » qa,qb is rather extreme. Instead, 
in the present scattering geometry, qc/cja, qc/qt will grad- 
ually decrease as the energy is increased. This will make 
the factors weighting the transverse and the longitudinal 
part to depend on energy. Due to the larger sensibil- 
ity to the transverse part in a real experiment this will 
lead to an apparent loss of spectral weight as the energy 
increases. 

In an ideal experiment where all components are de- 
tected with equal sensitivity one will not see a loss of 
spectral weight but rather a transfer of spectral weight 
from the longitudinal to the transverse part as the energy 
increases. Indeed S'^^-^{q,oj) satisfies the sum rule 

1 fOO 

^o" ^mT. d{nu;)S^ff{q, u) = 2(n - 27^)^1 

q 

(17) 

independently of what the distribution of domains or the 
scattering geometry is. 



D. Experimental Considerations 

On analyzing experimental data one should take into 
account that unpolarized INS experiments magnetic and 
nonmagnetic scattering can not be unambiguously identi- 
fied. In the insulating phase this problem can be reduced 
by fitting the data with a simple model like spin-wave the- 
ory. Due to the lack of a simple theory in the doped phase 
a similar procedure is not possible. In this case the usual 
experimental practice is to report as "magnetic scatter- 
ing" only those features which satisfy certain criteria like 
being reasonably sharp in momentum space and show a 
"magnetic like" temperature dependence and form factor 
dependence across different Brillouin zones. All the rest 
is assumed to be background of unknown origin. Below 
we analyze the distribution of weight in the insulator to 
get a hint of the distribution of spectral weight in the 
different channels and where to expect sharp and where 
broad features. 



III. HALF-FILLING 

A. Spectral weights in the Heisenberg Model 

We start by neglecting the shielding factors and esti- 
mate the spectral weights in the Heisenberg model within 
spin- wave theory (SWT). We will show below how to ap- 
ply these results when shielding factors are important 
corrections. Since we are interested in gross distributions 
of spectral weight for simplicity we consider a Heisen- 
berg model with nearest neighbor interactions and ne- 
glect other terms like four-ring exchange. 



1. Longitudinal part 

In order to write the longitudinal dynamical structure 
factor we introduce the following notations. For an AFM 
the operator can be written within SWT as: 

where A/i is the number operator for Holstcin-Primakoff 
bosons and Qafm = (■"', "■) taking the lattice constant 
a = 1. For the Fourier transform = J^i e^^'^ '^'S^ we 
have: 

5| - NSS{q - Qafm) - Mq-Q,,^M 

It is convenient to write N'q-Q^pM = {^q-Qafm) + 
^^q-QAFM with (A/q-Q^j^j,,) = NASS{q - Qafm) and 
define rnq^pj^j = m = S — AS as the reduced sublattice 
magnetization. Here 



N ^2 
q 
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is the reduction in the ordered spin moment due to zero- 
point quantum fluctuations and 7q = [cosq^ + cosgy]/2 
comes from the Fourier-transform of the exchanges in the 
square lattice. 

The one-magnon dispersion relation is given by: 

fiujq = 2JZc^l-ll (18) 

where J is the superexchange constant and Zc is 
a quantum renormalization of the one-magnon en- 
ergy, near-constant over the Brillouin zone in a first 
approximatioufiSiii and estimated as 1.18 to order 1/5^ 
in spin- wave theorjiiS or by series expansions.^"' 




(37c/2,jt/2) (rc,0) (3;:/2,x/2) (2j:,0) (rc,0) 



Wavevector {q ,q ) 

FIG. 1: Two-magnon scattering intensity [Eq. 1201 1 as a func- 
tion of energy and wavevector along symmetry directions in 
the Brillouin zone (solid lines in Fig.|21inset). Density of scat- 
tered points represents intensity. The lower boundary traces 
the one-magnon dispersion relation Tujjq and corresponds to 
events where one of the two magnons has zero energy. The 
upper bound, iZcJ, is reached when both magnons are on 
the antiferromagnetic zone boundary contour (dashed square 
in Fig. 121 inset) where they have maximum energy. 

With these definitions we can write the longitudinal 
dynamical structure factor as: 

SfiJq,u;) = g^yl[NmH(q-QAFM)5{hLo) (19) 

where the first term in the brackets is the Bragg elastic 
contribution of the Neel order and the second term is the 
inelastic contribution. 

The sublattice magnetization is well known to be accu- 
rately given by linear SWT^ therefore we do not expect 
significant changes on the elastic intensity if higher order 
corrections are included. In contrast how the inelastic 
part is distributed at low-energies (a few J's) is expected 
to be sensitive to such corrections. 

The inelastic part can be decomposed into a sum over 
inelastic processes out of which the dominant term cor- 
responds to two-magnon scattering events with intensity 




(371/2,31/2) (Ji.ji) (Ji,0) (3ji/2,ji/2) (2ji,0) (3i,0) 



Wavevector {a ,q ) 

X y 

FIG. 2: Energy integrated spectral weight in the two-magnon 
continuum (dashed line), S^mW) ~ J d{TiLo)S2M{q,i^), com- 
pared with the weight in the one-magnon peak (solid line), 
from Eq. We have used Z^ = 0.57 and Z2M = 0.67. 

given bjii& 

SriM,^) = Z^M^^ (20) 

X X! fiQi,Q2)Sihu; -hwq^ - hujq^)6{q + qi -q-a)- 
qi,<j2 

The two-magnons have opposite spin Sz — +1- and — 1 
such that the total spin is unchanged. The scattering 
cross-section f{qi, 92) depends on the wavevectors of the 
two magnons viei^ 

/(qi, 92) = sinh^ (6lqi -Og,_), tanh"^ 7,, 

(21) 

We have included an ad hoc intensity renormalization 
factor Z2M to be discussed below. The lineshapes in 
Eq. (|20|l can be evaluated by direct summation or Monte- 
Carlo methods and the result is illustrated by the dotted 
areas in Fig. Q 

Two-magnon events contribute a continuum band of 
scattering at energies above the one magnon dispersion 
TiLOq. In general the intensity decreases with increasing 
energy and it is strongest for wavevectors near the anti- 
ferromagnetic zone center (tt, tt). 

Fig-Elshows how the energy- integrated intensity varies 
in the Brillouin zone, generally following the same trend 
as the one-magnon intensity. The two-magnon signal 
cancels at the nuclear zone center (27r, 0) as expected 
from first-moment sum rules for an isotropic antiferro- 
magnet, also made apparent from Eq. 1)21(1 where the 
cross-section for creating two identical magnons cancels, 
i.e. for qi = q2 such that g = 0, /(<7i, gi) = 0. 

Fig.|21shows the wavevector-integrated intensity: it in- 
creases with increasing energy, linearly at low energies, 
then reaches a maximum at the one-magnon zone bound- 
ary and then decreases to zero at energies above twice 



the maximum magnon energy. In fact most of the two- 
magnon scatterig weight is at energies above the one- 
magnon cutoff 2ZcJ- 

Due to the broadness of the distribution the two- 
magnon spectral weight will be quite hard to detect in 
unpolarized neutron experiments, although we note that 
recent neutron experimentsii have observed such a high- 
energy continuum of excitations in the square-lattice 
spin-5/2 Heisenberg antiferromagnet Rb2MnF4 and the 
observed intensities were in agreement with neutron scat- 
tering by pairs of magnons as described by spin-wave 
theory. In cuprates, in addition, a large fraction of the 
spectral weight is at energies that are too high for neu- 
tron scattering. 

Integrating Eq. H20|l over energy and wavevector in 
a Brillouin zone gives Z2m^S{1 + AS) whereas 
the Bragg elastic scattering is (5 — A5)^ ff^Ms- The 
magnetization reduction can be evaluated numerically as 
AS* — 0.197. The corresponding value of Afg is shown in 
Table m Comparing with the Ising AFM, also shown in 
Table m one sees that the effect of the transverse fluctua- 
tions is to reduce the sublattice magnetization with a con- 
comitant transfer of spectral weight from the Bragg peak 
to the two-magnon continuum. Most of the longitudinal 
scattering is in this spread-out two-magnon continuum. 
The spin-only sum rule Mq = is strongly violated if 
intensity renormalization is neglected {Z2M = !)■ Notice 
that if transverse fluctuations are neglected (Ising limit) 
the sum rule is exactly satisfied. We can anticipate a 
similar result in the Hubbard model. 

Including higher-order terms in spin-wave theory is 
expected to produce: i) intensity-lowering of the two- 
magnon response^^ and ii) spread out at higher energies 
due to the contribution of higher multimagnon process^ 
For simplicity we assume the intensity renormalization of 
the two-magnon response to be momentum independent 
and quantify it with the constant Z2m- 

Since the Bragg intensity is given accurately by SWT 
one can get an upper bound for Z2M by neglecting effect 
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FIG. 3: Momentum integrated spectral function of one- 
[Eq. 11221 1 and two-magnon excitations [Eq. I20ll 1. We have 
used Zd = 0.57 and Z-^m = 0.67. Integrated areas are 
(1 2/^S)ZiS and Z2m^S{1 + AS), respectively. 
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ii) and requiring that the sum rule Eq. Q is fulfilled by 
Bragg and two-magnon processes: Z2M ^ 0.67. Pertur- 
bative computations^ suggest a value not far from that 
suggesting that the weight in four magnon and higher 
multimagnon processes is not high, in agreement with 
numerical data^ In the following to improve the val- 
ues of the two-magnon intensity we tentatively adopt 
Z2M = 0.67 given the lack of better estimates (Table QJ. 

It is interesting to remark that the small weight at four 
magnon and higher multimagnon processes is in striking 
contrast to the response relevant for infrared experiments 
which show large weights insteadi^°i^-^i^^i^'^i^^'^^ From a 
theoretical point of view, whereas magnon-magnon inter- 
actions have dramatic consequences in the shape of the 
spectrum of two-magnon excitation relevant for Raman 
and IR data^-'^'^^ the effect on the two-magnon line shape 
relevant for magnetic neutron scattering appears much 
modest This different role of interactions can be traced 
back to the fact that optical data probes two magnons in 
different sublattices whereas magnetic neutron scatter- 
ing probes two magnons on the same sublattice and the 
interaction is dominant in the first case. The correspond- 
ing Green functions at the RPA level for both cases are 
reported in Ref. 22. In optical data magnon-magnon in- 
teractions are responsible for a shift of the main feature 
from ~ 4J to 37,21,22^26 rpj^jg effect was obtained in 
Ref. 22 using a high-energy approximation. One can eas- 
ily see from the RPA equations for the equal-sublattice 
Green functionS^ that at the same level of approximation 
the same effect does not show up in the two-magnon neu- 
tron scattering line shape. We conclude that interaction 
effects should be important for the intensity renormaliza- 
tion but they play a very different role here than the one 
played in optical data. 



TABLE I: Sublattice magnetization mg^j^j,^ and dynamical 
structure factor spectral weights in units of fjig . We show ex- 
act values for the Ising AFM and SWT values for the Heisen- 
berg AFM in the longitudinal channel (2) and in one trans- 
verse channel (x). SWT are the values neglecting all renoma- 
lization factors {Z = 1). We also include the SWT intensity 
renormalized by the indicated Z values. The transverse con- 
tribution is split in the contribution from the magnetic zone 
and the nuclear zone as Mfigz+Mfjgz- Since we are consid- 
ering the Heisenberg model {D = 0) a "1" in the last column 
implies that the sum rule is exactly satisfied within the model. 
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Inelastic 
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0.5-h0.5 
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SWT 
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0.3034 


0.368 


0.941" 


1.309 
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1.068-^0.325* 
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Z2M = 0.67 
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0.3034 


0.368 
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Za = 0.57" 


X 







0.609+0.185' 


0.794' 



"Two-magnon contribution. 
'One-magnoncmitribution. 
After Refs. Iljld 
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2. Transverse part 

The transverse dynamical structure factor is domi- 
nated by one-magnon scattering events with intensity 
given bjii^ 



2 2 r, 

9 Aifi Zd 



S 1- 



7q 



1-7^ 



5{nuj - hujq), (22) 



Zd is an intensity-lowering renormalization factor of the 
one-magnon cross-section due to zero-point fluctuations 
and magnon-magnon interactions, both neglected at first 
order in spin-wave theory [Zd — Z^Zc, where Z^ is the 
renormalization of the transverse magnetic susceptibility, 
X± = Z^ig^iB)VSJ). 

The wavevector dependence of the one-magnon inten- 
sity is shown in Fig. 12 Although wave-vectors {2tt, 0) = 
(0, 0) and (tt, tt) are related by symmetry the spectral 
weight goes to zero at wave vector {2tt, 0) and diverges 
at (tt, tt). This difference in intensity is reflected also in 
doped phases. 

Integrating Eq. (|22|) over energy and wavevector in a 
Brillouin zone gives the one-magnon intensity for one 
transverse direction as (1 + 2AS)ZdSg^iJ,^/2. Experi- 
mental works often restrict to the magnetic Brillouin zone 
around (tt, tt) (dashed line in the inset of Fig. 12) hereafter 
"the magnetic zone" and neglect the small weight on the 
magnetic Brillouin zone around (0, 0) hereafter "the nu- 
clear zone" . For comparison we split the one-magnon 
spectral weight in the two zones Mg = M'fj^^ + ^^nbz- 
One finds the values reported in Tabled 

In this case as in the longitudinal channel the sum rule 
is overestimated if one neglects the intensity renormaliza- 
tion. Estimates of the renormalization from higher-order 
spin-wave theory giveiSii^ Zd = 0.57 at order 1/5^, or 
0.61(4) by series expansions^ Now the sum rule is un- 
derestimated. The lacking weight is expected to lay in 
three-magnon and higher multimagnon processes. 

An alternative way to estimate Zd is to enforce Eq. Q 
at large S and extrapolate to S" = 1/2. With the choice 
Zd = 1- AS/S = 0.606, the total sum rule S{S + 1) 
[c.f. Eq. (Q] is exhausted by elastic {S — AS*)^, one- 
magnon (1 -I- 2AS){S — AS), and two-magnon scattering 
AS'(H-AS') (without renormalization), in units oig^jj^g. 
Such a renormalization of the one-magnon intensity is 
very close to the results quoted above. Strictly speaking, 
this argument applies only for S" > 1/2 since the sum rule 
for each channel Eq. ©, which is exclusive of spin 1/2 
systems, is still strongly violated. One the other hand the 
rapid convergence of the l/S expansions often produces 
accurate results for S — 1/2 systems as seem to be the 
case for Zd- 



B. Shielding factor in the Hubbard Model 

We now turn to the more realistic Hubbard model. In 
this section we estimate the shielding factor for n = 1, 



within the 2d single-band Hubbard model Eq. Q , to gain 
some insight on its impact on the sum rule Eq. © . Here 
we take for simplicity t' = and consider the effect of 
varying U /t. Parameters more specific for the cuprates 
will be considered in Sec. IIVI 

The shielding factor is given by the single occupancy 
probability n — 2D. For a non interacting system, i.e. 
[/ = 0, the double occupancy is just (n/2)^ so one gets 
n — 2D — 1/2. Indeed half of the time an atom is singly 
occupied (up or down) and half of the time it is either 
empty or doubly occupied and hence does not produce 
scattering. 

For a paramagnetic state in the Hartree-Fock approx- 
imation (HE) one gets the same result independently of 
U /t (Fig. ^ since correlations are neglected. Correla- 
tions can be introduced by treating the paramagnet in 
the Gutzwiller approximation (GA)^ which leads to a 
reduction in double occupancy as 



4 V [/, 



becoming zero at the Brinkman-Rice transition point 
Uc/t = 128/7r2 ~ 12.97,^2iSa The corresponding value 
of n — 2D is plotted also in Fig. 01 The Brinkman-Rice 
point, however, is never reached since for infinitesimal U 
the paramagnetic state is unstable towards a spin-density 
wave (SDW) which can also be treated in the OA and 
gives the value of n — 2D shown in Fig. ^ 

Finally for large U /t one can use a canonical trans- 
formation to map the Hubbard model to a Heisenberg 
modcl;i^*2S It is important to realize that the energies of 
the low-energy excitations of both models, the "physical" 
Hubbard model and the low-energy "effective" model co- 
incide to leading order but the correlations functions in 
general do not coincide. In order to get "physical" cor- 
relation functions one needs to use the inverse canonical 
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FIG. 4: (Color online) Single occupation probability (shield- 
ing factor) for a half-filled Hubbard model in different ap- 
proximations. The exact results in a 4 x 4 site cluster are 
computed with the double occupancies given in Ref. W\ The 
apparent lack of extrapolation of the exact results to the non- 
interacting limit at f/ = is due to a finite size effect. 
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transformation to transform back the "effective" ground 
state wave function to a physical wave function. Indeed 
within the Heisenberg model the double occupancy is 
zero but this does not mean that the double occupancy 
is zero in the "physical" model. This is obviously very 
important on evaluating the right hand side of Eq. JH)). 
Of course if we are interested on evaluating the sum 
rule within the Heisenberg model it is legitimate to take 
D = 0, as done above, but if we want to compare with 
experiments (or with the "physical" model) one needs to 
compute the "physical" double occupancy. Fortunately 
this is very easy in the present case because we can use a 
trick to avoid the back transformation. We use Hellman- 
Feynman theorem to write the double occupancy in the 
Hubbard model as: 

N dU 

with E the ground state energy. For the latter we use 
the fact that for large U /t it coincides with the energy 
of the Heisenberg model where very accurate numerical 
estimates exist: 



E 



it' 



and the best estimate for a is a = 0.6696ii^ Here we have 
substituted the superexchange constant by its definition 
J = At'^/U. Since the Heisenberg model has only one 
parameter it is clear that a does not depend on J and 
we can perform the derivative to obtain: 



D 



At' 



The corresponding value of the single occupancy prob- 
ability n — 2D is plotted in Fig. ^ and referred to as 
the H/U expansion". This value is asymptotically exact 
(within numerically accuracj*i^) in the large U jt limit 
and therefore takes into account all fluctuations effects. 
Clearly the SDW treated within the GA approximations 
gives a fairly good approximation for n — ID at large V 
and interpolates smoothly to the exact result at C7 = so 
we expect it to be quite accurate in all the range of U jt 
as can be seen also by comparing with the exact results 
in a 4 X 4 cluster (after Ref.ll^). 

We see that the sum rule for one diagonal component of 
the dynamical structure factor Eq. (O changes smoothly 
from in the non- interacting case to /z^ in the limit 

of [/ = oo. It is reduced by - 11% (Table HJl with 
respect to the full moment value for [//< ~ 8, as relevant 
for cupratesi^Si^ 



Distribution of spectral weight in the Hubbard 
Model 



In the previous section we have evaluated the shielding 
factors that appear due to the finiteness of V and reduce 




FIG. 5: (Color online) Dispersion relation of the low-energy 
transverse excitations. We show the experimental result for 
La2Cu04 after Ref. ^ and the GA-I-RPA resuh for Vji = 8 
and t = 335meV. The inset shows that for U/t = 10 the fit is 
noticeably worse. 



the total spectral weight compared to spin only models. 
Another effect, which we study here, is that the total 
weight is split into a low-energy part at energies of order 
J and a high-energy part at energies of order U. Present 
neutron scattering facilities can measure the spectra up 
to energies of the order of a few tenth of eV and therefore 
only the first part is detected. In addition we discuss 
how the dynamical structure factor is modified for finite 
U respect to the SWT result. 

To estimate the dynamical structure factor we use 
the time dependent Gutzwiller approximation of Refs. 0, 
[also called GA plus random-phase-approximation 
(GA-hRPA)] applied to a SDW state. We start in the 
next section by showing how measurements of the dis- 
persion relation in the insulator can be used to estimate 
U and t. 



1. Transverse part 

In the transverse channel MNS experiments reveal 
the spin-wave excitations of the AFMi^ Interestingly, 
whereas spin-wave theory on the Heisenberg model 
Eq. (|18|l predicts a flat dispersion between (tt, 0) and 
(7r/2, 7r/2) a substantial dispersion has been measured 
(c.f. Fig. EJ. It has been argued that in cuprates cor- 
rections to the Heisenberg model arising as higher orders 
in a t/U expansion are relevant iSSi^SiSSiSi^ The most 
important of such corrections is a term which cyclically 
exchanges four spins on a plaquette. A sizable value for 
this term has been revealed by analyzing phonon-assisted 
multimagnon infrared absorptions^ and the dispersion 
relatioii^ shown in Fig. In particular the dispersion 
between (tt, 0) and (7r/2, 7r/2) is mainly due to this term. 
Since the dispersion has its origin in the finiteness of t/U 
it should show up in the transverse excitations of the 
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Hubbard model. The computation done in GA+RPA is 
also shown in Fig. [S] One obtains a very good fit of 
the dispersion and this provides an accurate way to esti- 
mate the strength of the repulsion. We find U /t = 8 in 
good agreement with other estimatesSA. whereas U /t = 10 
gives a too flat dispersion between (tt, 0) and (7r/2,7r/2) 
(inset). An equally good fit as the one show in the main 
panel can be achieved with U/t = 8, t'/t = —0.2 and 
t = 353.7meV. The value of t'/t plays an important role 
in the doped phasa^ii^ and is close to a first principle 
estimate 1^ 

In Fig. we show the momentum integrated spectral 
function. We see that for large U/t it approaches the 
SWT form (Fig.O. For smaller U/t the maximum is not 
at the upper edge due to the modified dispersion relation. 
A small portion of the spectral weight is at high-energy 
due to spin-fiip transitions from the lower to the upper 
Hubbard band as shown in the inset. We will show that 
this effect is much larger in the longitudinal channel. 

Integrating the spectral intensity one gets the values of 
Mq shown in Table Ull GA-I-RPA interpolates between 
the extreme limits U/t = (where the sum rule is exactly 
obeyed) and U —>■ oo where one recovers the results of 
linear spin-wave theory. As a consequence the sum rule 
in the transverse channel is increasingly more violated 
as U/t increases with Mq reaching values similar to the 
unrenormalized SWT results (Table HJ for large U. 

In Fig. [3 we show the experimental intensities as a 
function of wave vector together with the GA-I-RPA re- 
sult. In order to fit the experimental results we in- 
troduced an intensity renormalization Z^^^* analogous 
to the SWT intensity renormalization factor. With 
= 0.65 the GA-I-RPA intensities are essentially 
equivalent to the intensities reported in Fig. 3b of Ref . 1301 
which have been renormalized by Z^^^ = 0.51 respect to 
LSWT. 

This increase in the value of Z^ reflects the fact that 
cuprates are not in the strict Heisenberg limit U —>■ oo. 




(7i/2,ji/2) (Ji.jt) (Ji,0) (7i/2,jt/2) (0,0) 

FIG. 7: (Color online) Intensity in the spin-wave excitations 
as a function of momentum for one transverse spin-channel. 
We show the experimental result for La2Cu04 after Ref. IscI 
and the GA-I-RPA result for U/t = 8 and t = 335meV. 
The GA-I-RPA intensity has been renormalized by a factor 
Zd^^* ~ 0.65 and coincides with the spin-wave result (shown 
in Fig. 3b of Ref.iSfl) renormalized by Z^""^ = 0.51 in Eq. ll^ . 



In that limit the Zd for the RPA should coincide with the 
value of LSWT given the equivalence of the two approxi- 
mations. In the extreme case of a non interacting system 
RPA gives exact intensities and therefore Z^^'^ = 1. As 
U is decreased from infinity the RPA becomes gradually 
more accurate and therefore Z^ should increase respect 
to the LSWT value as we indeed find. 

We caution that the fact that a system is not in the 
strict Heisenberg limit does not mean by itself that it 
can not be described by the Heisenberg model. Correc- 
tions to the intensity can be introduced as explained in 
Sec. nil Dl We estimate this approach is accurate up to 
U ^ 15t where the dispersion is still Heisenberg like. For 
smaller U corrections in the dispersion relation due to 
four spin exchange and other non Heisenberg processes 
become obvious (c.f. Fig.lSj). 




FIG. 6: 5*^^(0;) evaluated for the Hubbard model within 
GA-f- RPA for U/t = 8 (solid line) and U/t = 10 (dashed line) 
with = 1 for an 80 X 80 sites system. The inset shows the 
high energy contribution of S^^{u!) at energies u ~ (7. 



2. Longitudinal part 

The longitudinal structure factor has a Bragg part and 
an inelastic part [c.f. Eq. H12|l ]. For the SDW only one 
magnetic reciprocal lattice vector contributes to the sum 
in Eq. H12() namely Qafm = (i", ti") (we take the lattice 
constant a ~ 1). The Bragg weight is determined by the 
sublattice magnetization rnq^pj^j . 

At lowest order RPA introduces longitudinal fluctua- 
tions in S^^ but does not correct one-body expectation 
values like the sublattice magnetization. (On the con- 
trary it introduces a correction in two-body expectation 
values^ like D). Therefore in the following we consider 
the sublattice magnetization at the GA level together 
with the effect of the longitudinal fluctuations at high en- 
ergies in S^^ . The sublattice magnetization can in prin- 
ciple be corrected by including transverse fluctuations. 
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We will discuss this below. 

At the GA we obtain the values of the sublattice mag- 
netization shown in Table HTl These values lack the trans- 
verse fluctuation corrections and hence converge to the 
classical value, 0.5 in the large U limit. 

To help intuition it is useful to distinguish between 
the permanent moment m — — 7iij^)|/2, which de- 
termines the Bragg weight, and an "intrinsic" moment 
which determines the shielding factor and which we de- 
fine as TO* = V(("»T - "a)^)/2 = A/(n-2D)/2. (No- 
tice that with this definition a free fermion system has 
an intrinsic moment which we denote as "trivial"). The 
paramagnet with small D mentioned in Sec. lIIIBl can be 
considered as a system with well formed (non-trivial) in- 
trinsic moments that are purely dynamical. Indeed the 
shielding factor is close to 1 but there are no Bragg peaks. 
In this regard the GA is more flexible than HF which is 
not able to produce non-trivial intrinsic moments that 
are not permanent. 

Permanent moments are due to the breaking of spin- 
rotational symmetry. The GA for the SDW smoothly in- 
terpolates between the itinerant limit at small U and the 
localized limit at large U with a permanent moment that 
increases and reaches 1/2 for f/ — > oo. Contrary to the 
paramagnetic case, quoted above, the difference between 
permanent and the intrinsic moment in the GA tends to 
be very small at large U. Thus despite the ability of the 
GA to distinguish the two kinds of moments this ability 
is not effective at large U. In this case the permanent 
moment in the GA (c.f. Table ITl|l is reduced with respect 
to the fully polarized value mainly due to covalency ef- 
fects. On top of this the moment will be reduced due to 
transverse fluctuations as discussed below. 

The covalency reduced Bragg peak does not exhaust 
the sum rule as can be seen by comparing the "Elas- 
tic" column with the last column in Table ^ Lon- 
gitudinal fluctuations, captured by the GA-fRPA ap- 
proach, produce weight at energies of order U as shown in 
Fig. |S| This is not related to two-magnon processes (as 
in Sec. IIII A l|l but to fluctuations in the length of the 



TABLE IL Spectral weights in units of fj,%. We show exact 
values for the Hubbard in the longitudinal channel (z) and in 
one transverse channel (a;). The inelastic weight is separated 
in the low-energy part (Low) at energy ~ J and the high- 
energy part Mo, (7 around energy ~ U. We also show the 
sublattice magnetization mq^^j^j and the shielding factor n — 
2D, both computed in the GA (without RPA correction). 
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FIG. 8; (Color online) Inelastic part of S^^(ijj) in the 
GA-I-RPA for different values of U/t in an 8 x 8 sites sys- 
tem. 



intrinsic moment. The weight in the high-energy part 
Mo,{7 is also reported in Table Hll Adding the Bragg and 
the Mq jj contribution one obtains the values of Mq re- 
ported in Table ITU Comparing with the last column we 
see that at this level of approximation the sum rule is 
exhausted (actually slightly overshot) by the Bragg con- 
tribution without transverse fluctuations reductions plus 
the contribution at energy of order U . 

The small overshot in the sum rule is due to the fact 
that the last column is computed at the GA level without 
RPA corrections, whereas Mq has RPA corrections to D. 
(An analogous method was used in Ref. to compute 
fluctuation corrections to D). The Mq column can be 
considered as an improved computation of the shielding 
factor. We see that the difference is very small and for 
practical purposes we can use the GA result. 

From the above results it is clear that neglecting trans- 
verse fluctuations the sum rule is well satisfled. This is 
analogous to the behavior found for the spin only model 
in the Ising limit (Table HJ. 

At this point we have weight at the Bragg peak and 
weight at energy of order U . How these weights get af- 
fected by transverse fluctuations? On top of the cova- 
lency reduction transverse fluctuations reduce the perma- 
nent moment (keeping the intrinsic moment constant). 
We have seen above (Sec. IIII A l|l that the reduction 
of the moment due to transverse fluctuations produces 
a concomitant large reduction of elastic spectral weight 
which is transfered to multimagnon processes in the lon- 
gitudinal channel at energy ^ J. A direct computation 
in the Hubbard model is difficult. In the next section we 
show how to estimate this effect in the Hubbard model 
using the SWT results. 

We expect that this low-energy rearrangement of spec- 
tral weights will neither affect Mq^ nor the shielding 
factor computed above, which, as is clear from Fig. ^ 
is quite accurate at the present level of approximation. 
Therefore we expect the present values of (n — 2D)fi^ — 
Mqjj and Mq^ (without transverse ffuctuations) to be 
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FIG. 9: (Color online) Momentum integrated longitudinal 
spectral function for the Hubbard model with U — 8t, t — 
335meV and ZcJ = 153meV. We show the contribution of 
two-magnon excitations computed in SWT [Eq. 12311 1 with an 
effective shielding factor of 0.74 (Table El and Z2M = 0.67. 
We also show the contribution due to scattering across the 
Hubbard bands computed in GA+RPA. 



accurate estimates of the total spectral weight at low (or- 
der J) and high (order U) energies. 



D. Effective shielding factor 

In this section we would like to show how to use the 
spin wave theory results together with the results of the 
previous section to obtain a better distribution of the 
spectral weight. 

According to the discussion of Sec lIII B1 to obtain the 
"physical" dynamical structure at low energies from the 
Heisenberg model response, one should apply the canon- 
ical transformation back from the effective model to the 
Hubbard model. Since the processes involved in a t/U 
expansion have a short spatial range and involve inter- 
mediate states at high energies, we do not expect this 
procedure to lead to a momentum or low-energy depen- 
dent correction. Thus the "physical" dynamical correla- 
tion functions at low energy can be obtained from the 
requirement that the sum rule is satisfied as 

5""(g,a;) - (n -2D- MoV//i|)^J^e.(9, c^) (23) 

where S'^^{q, lo) is the dynamical correlation function of 
the Heisenberg model and M^^ is the weight transfered 
to high energies due to the finiteness of U computed in 
Sec. nil CI Alternatively a model with higher order cor- 
rections (ring exchange, etc.) can be taken to compute 
the response on the right hand side of Eq. H23|l . 

We can consider {n -2D- M^jj/fj^g) in Eq. as an 
effective shielding factor for low-energy spectral weight. 
One can check from the results of TablelTllthat this works 
for the Bragg intensities neglecting transverse fluctua- 
tions. Indeed we have 



where the left hand side is the Bragg weight in the GA 
for the Hubbard model without RPA correction and the 
right hand side is the effective shielding factor times the 
Bragg weight in the Ising limit of the spin only model 

Due to the separation of energy scales D and Mq ^ will 
be insensitive to transverse fluctuations. The effect of 
the latter will be to transfer weight from the Bragg peak 
to the low-energy continua without affecting the total 
low-energy spectral weight. Thus Eq. tell us that 
the total low energy spectral weight (Bragg plus multi- 
magnon contribution) is accurately given by the mean 
field value 4mo Mr- 

As a corollary we see that a good estimate of the longi- 
tudinal effective shielding factor is given by the moment 
in the GA without the need of an RPA computation, 
i.e. the low-energy part of the longitudinal dynamical 
structure factor gets rescaled by the mean-field value of 
4:171% . This is consistent with the findings of Hirsch 
and Tang from numerical data™ if one takes into ac- 
count that the HP approximation, used by Hirsch and 
Tang, gives very similar values of toq^^m as the GA in 
the range of U analyzed. 

In Fig. El we show the two-magnon contribution com- 
puted in the Heisenberg model and translated into a Hub- 
bard response according to Eq. H23|l . together with the 
longitudinal response at energy ~ U computed directly in 
the Hubbard model. Additionally there is a Bragg peak 
(not shown) with weight 16mQ^^^j/i^(S' — AS)'^ where 
^Qafm is computed at mean-field level and 4(5* — AS*)^ 
takes into account the effect of transverse fluctuations. 
This is our best estimate for the distribution of spectral 
weight in the longitudinal response of the insulator. 

Notice that the transition across the Hubbard bands 
appears at 2eV. One should remember that on the map- 
ping from the three-band Hubbard model to the one- 
band Hubbard model2& the transition across the Hub- 
bard bands of the latter represent the charge transfer 
transitions of the former, which according to optical data 
should occur close to ^ 2eV.^S Indeed in the ionic limit, 
the lower Hubbard band of the one-band model corre- 
sponds to O states and the upper Hubbard band to 
Cu dio. 



E. Estimation of the weights in La2Cu04 and 
Cu(DCOO)2.4D20 (CFTD) 

In La2Cu04 the dynamical structure factor can be 
fitted^ in a wide range of momentum and energy by 
the SWT expressions Eq. H22|) with an intensity-lowering 
renormalization factor = 0.51 ± 0.13. We can use 

this result to estimate the different contributions to the 
moments for the one-magnon processes in La2Cu04: 

M^B2 = 0.16 ±0.04^*1, 



4m: 



Qafm 



2D)f,l 



M, 



o.u 



(24) 



0.54±0.13Ai| 
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and 

= 0.71 ± 0.18^1%. 

The last value should be interpreted as the weight for one 
transverse direction in the spin-wave like excitations of 
La2Cu04. This is not the total weight since part of the 
weight will be in higher multimagnon processes. 

The experimental determination of the intensity low- 
ering factor, Z^^^ neglected shielding factors. Therefore 
the experimentally determined quantity can be put as: 
^e.p ^ ^eff ^ _ _ MS^u/t^B)Zd. Using an ef- 
fective shielding factor for U/t = 8 as appropriate for 
La2Cu04, {n-2D - M^jj/^jlb) = 0.88 and Zd = 0.57 

(SWT) one obtains Z'^^^ — 0.50 in agreement with the 
experiment (inclusion of t' /t — —0.2 does not change this 
number appreciably). 

Subtracting the observed one-magnon weight from the 
total expected transverse weight of 0.88/is would then 
leave a total of 0.17 ± 0.18/i^ in three, five and higher 
multi-magnon processes which is within the error bars 
as illustrated in Fig. ^| Notice that the detected one- 
magnon spectral weight is significantly lower than 1/xb, 
the value expected from the sum-rule neglecting the 
shielding factors and quantum corrections. 

For Cu(DCOO)2.4D20 (CFTD) R0nnow et al^ also 
find Z^^'P = 0.51 ± 0.04 within LSWT. For this com- 
pound the ratio of U /t is much larger than in the cuprates 
and the shielding factor correction should therefore be 
close to 1. Since the error bars are smaller than in the 
cuprates this experiment shows, not surprisingly, that 
multimagnon processes are needed to satisfy the sum 
rule (see Fig. llOfl. The similarity for Z^'^^ found for the 
two compounds is to some extent unexpected since there 
should be a difference due to the shielding factors, how- 
ever, the expected difference is within the current exper- 
imental errors. Notice also that g may also differ in the 
two compounds. 

IV. AWAY FROM HALF-FILLING 
A. Shielding factors 

Taking g ^ 2 and defining the doping as a: = 1 — n, 
the shielding factor for general doping changes from (1 — 
a;^)/2 for U = to 1 — \x\ for U = oo, the latter form 
being used often in experimental works. In general we 
expect a dependence on U/t qualitatively similar to the 
one shown at a; = (c.f. Fig. 0. 

In order to proceed in the doped phase we need a model 
for the ground state. In this work we are interested in 
overall distributions of weights and we expect this to be 
to a large extent insensitive to the details of the ground 
state if correlations are taken reasonably well into ac- 
count. This can be seen already at x — (c.f. Fig. ^ 
where we see that the paramagnet and the SDW within 
the GA give similar shielding factors although the nature 



of the ground state is completely different. Detailed dis- 
tributions of spectral weight as measured in Ref. Q and 
computed in Ref. |^ will, of course, depend strongly on 
the ground state and so can be quite helpful to deter- 
mine it. Since those issues are beyond our present scope 
we restrict the description of the mean-field states to a 
minimum. 

For the specific system La2_2;Xa;Cu04 (X=Sr,Ba) we 
use a ground state consisting of stripes as suggested by 
experiment Metallic stripes parallel to the CuO bond 
can be obtained within the GA for the Hubbard model 
where a next-nearest neighbor hopping t' /t = —0.2 has 
to be implemented in order to have one doped hole per 
every second unit cell along the stripe^^ (so called "half- 
filled" stripes) in agreement with experiments. Note that 
a similar value for t' /t is predicted by first principle 
computations'^^ for La2_a;Sr2;Cu04. 

It should be mentioned that the results at half-filling 
reported in the previous sections are rather insensitive 
to the next-nearest neighbor hopping. For U/t — 8 and 
t' /t = —0.2 we obtain an analogous dispersion than that 
reported in Fig. [S] with the only difference that the 
value for the nearest neighbor hopping has to be scaled 
to t = 353.7meV. This latter parameter set has been 
used^- in order to explain a recent MNS experiment on 
Ba-codoped LSCO.^ Overall spectral weights with this 
extended parameter set are practically the same than 
those reported in Sec. IIII (TTI with t' = 0. 

In Fig. El we show the shielding factor obtained aver- 
aging (n — 2D) in the mean-field solutions as a function 
of doping. We see that the shielding factor has a lin- 
ear behavior similar to the one for the U/t = oo case. 
As expected we find that the shielding factor is similar 
for other low-energy textures (not shown). For doping 
x = 1/8 we obtain n — 2D = 0.79 so there is an extra 
10% reduction with respect to the undoped case. 



B. Spectral weights 

1. Longitudinal Part 

Because the stripe solutions are magnetic the longitu- 
dinal structure factor has a Bragg part and an inelastic 
part [c.f. Eq. fT^ ]. 

Since the stripe solutions are metallic, low-energy 
particle-hole excitations are allowed. Therefore con- 
trary to the result in the insulator (Sec. IIII ("2(1 . al- 
ready at RPA level one finds weight at magnetic energies 
(^ 0.3eV) as shown in Fig. ^2 This weight however is 
broadly distributed in momentum space and will most 
likely pass unnoticed in unpolarized INS as discussed in 
Sec. Ill Dl Just as shown in Fig. for the insulator a 
more elaborate computation, taking into account longi- 
tudinal fluctuations, will show in addition weight at the 
two-magnon excitations. Since the inelastic longitudi- 
nal component due to multimagnon scattering is already 
quite featureless in the insulator f Sec. IIII A ll and FlII G 2|) 
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FIG. 10: (Color online) Zeroth moment and shielding fac- 
tor as a function of doping. We show the shielding factor 
for U/t = 8 and t' /t = -0.2 for the SDW at a; = and 
bond-centered metallic stripes solutions as a function of dop- 
ing (red line) and for U/t — oo (dotted blue line). For the ze- 
roth moment MBZ means that the integration was restricted 
to the magnetic Brillouin zone. La2Cu04^° and CFTD^ 
label the moments estimated in insulators (Sec. IIIIEII and 
La2-a;Baa;CuO^ labels the experimental value of Mo in the 
magnetic Brillouin zone, corrected by the polarization factor. 
In the last case the error bars reflect the two extreme possibil- 
ities for the polarization factor 2/3 and 1/2 (c.f. Sec.^. The 
relative error should be larger than at zero doping where the 
error bars have the usual sense. Finally we show the zeroth 
moment of the GA-I-RPA transverse spectra at j: = 1/8 in the 
whole Brillouin zone (-I-) and in the MBZ x. The two lower 
X include a Z^^^^ =0.65 renormalization factor and for the 
lower one the energy integration was restricted to energies 
smaller than 0.23eV. 



it is reasonably to expect that it will be even more fea- 
tureless in the doped phase. We conclude that within 
standard protocols (Sec. Ill Dp to first approximation all 
the inelastic longitudinal spectral weight will be assigned 
to the background. 

For the specific case of a striped ground state we can 
also estimate the elastic spectral weight in the longitudi- 
nal channel. For doping a; = 1/8 and an array of d = 4 
SC stripes we find elastic peaks at Qi = (1/2 ± e, l/2)27r 
and Q2 = (1/2 ± 2e, I/2)27r with mq^ = 0.215 and very 
small weight, TOq, in the higher harmonics mq^ = 0.018. 

It was shown in the insulator that if one neglects trans- 
verse fluctuations the sum rule in the longitudinal chan- 
nel is exhausted with a good degree of approximation by 
the Bragg weight plus the inelastic high energy weight. 
Is there an analogous behavior in the doped phase? The 
elastic weight neglecting transverse fluctuations is 0.37/x^ 
and the inelastic longitudinal spectral weight is 0.59/i^. 
This includes low-energy particle-hole excitations and 
weight due to transitions across the Hubbard bands. The 
total weight in the longitudinal channel is 0.96/i^ to be 
compared with a shielding factor of 0.79. In this case we 
flnd an overstrike of the sum rule which is more sever 
than in the insulator and which is even larger than the 
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FIG. 11: (Color online) Longitudinal and transverse com- 
ponents of the spin autocorrelation function for a; = 1/8 
in a 16 X 4 site cluster for U/t — 8 and t'/t = —0.2 and 
t = 353.7meV. Only the inelastic part is shown and we ne- 
glect the intensity renormalization {Z — 1). The ground state 
consist of 4 SC stripe running along the short dimension. 



maximum allowed value 1 — x ~ 0.875. This is because 
RPA becomes less accurate in the metallic phase due to 
the presence of small energy denominators. The violation 
however is not as severe as for the transverse channels as 
reported below. 

The elastic weight neglecting transverse fluctuations is 
roughly half of the one found in the insulator (0.74/iB)- 
Just as in the insulator a large fraction of this weight 
(at least 60%) will be transfered to multimagnon exci- 
tations if transverse fluctuations are taken into account. 
It is conceivable that all the weight in the Bragg peak 
is transfered in which case the system becomes quantum 
critical or quantum disordered^S 



2. Transverse Part 

The transverse component of the dynamical structure 
factor shows sharp features due to the propagating spin- 
wave like modes of the stripes. In addition, since the 
stripes are metallic, a fraction of the spectral weight is 
in a broad particle hole continuum. The latter features 
however are mainly located in the nuclear Brillouin zone 
(NBZ).3J^ 

In the magnetic Brillouin zone (MBZ) we flnd mainly 
sharp propagating collective modes which are good can- 
didates to be easily detected. These collective modes 
are analogous to the spin-wave modes of the insulator 
and therefore we expect that their spectral weight will be 
overestimated in RPA and should be corrected by a quan- 
tum renormalization as in the insulator. may be 
smaller than in the insulator due to stronger quantum 
fluctuations. On the other hand for a heavily doped sys- 
tem one should recover the noninteracting value Z^ = 1- 
Since doping is small we tentatively take the same value 
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FIG. 12: (Color online) The upper panel shows the theoretical 
and the experimental result for S'^^^ (uj) using U/t — 8, t' /t — 
-0.2, t = 353.7eV and (1 - ql))dom = (1 - qy))dom = 1/2, in 
a 96 X 96 sites lattice, together with the experimental data for 
La2-a;Baa;Cu04 X = 0.125* . The lower panel shows the par- 
tial sum rule weight. M^^^iuj) = j^J^g I-o d(huj')S^^ ' {q,u') 



as in the insulator Z^~^ — 0.65. 

The momentum distribution of the transverse spectral 
weight has been shown in Ref. 31. In Fig.^|we show the 
transverse integrated spectral weight compared with the 
experimental results. At zero doping the spectral weight 
coincides with the theoretical value (not shown for clar- 
ity) as discussed in Sec. lIIII At finite doping the weight in 
the MBZ decreases roughly as the shielding factor as can 
be seen by comparing the points labeled La2Cu04(MBZ) 
and RPAxZy=**(MBZ). Comparing Fig. [§1 and Fig. [TTl 
we see that whereas in the insulator almost all of the spec- 
tral weight is at magnetic energies in the doped phase a 
large fraction of the spectral weight is at intermediate 
energies not accessibly to MNS. In other words doping 
induces a transfer of spectral weight from low energies to 
intermediate energies in this channel. Interestingly opti- 
cal spectra shows that doping also generates structure at 
the same energies in the charge channel^^iiMi. 

The lower "x" at x = 1/8 in Fig. [101 is the RPA mo- 
ment in the MBZ renormalized by Z^^^^ and restricting 
the integral only up to the highest energy measured in 
Ref. Q, u> < 0.23eV. We see that the obtained value is 
roughly 50% higher than the experimental weight. 

In order to see at what energy the disagreement arises 
we show in Fig. E| the experimental S'^-^^ compared 
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FIG. 13: (Color online) Experimental results for S'^-^^{uj) 
for X = (Ref. 53), La2-^Ba^Cu04 x = 0.125 (Ref. [J 
and La2-a;Sr^Cu04 with x = 0.14 (Ref. 14^1) and x = 0.16 
(Ref. li^ . We used the following conversions from the quan- 
tity reported in the quoted reference: S"-^-^ = S (Ref. 0) 
= |x(Ref.|4i = fx (Ref.Ei. 

with the theoretical one using 1/2 as polarization factors 
[Eq. (|16|l ] and neglecting the longitudinal weight. The 
lower panel shows a comparison of the partial sum rule 
weight with the two limiting values of the polarization 
factor. The theoretical computation shows a strong res- 
onance at w = 58meV which in the experiment appears 
with a reduced spectral weight and in part distributed 
on a broader energy range. (The resonance appears at 
higher energy if one looks only at the Qafm responseSi). 
The overall weight integrated up to energies immediately 
above the resonance is in agreement with the experimen- 
tal one. As the energy increases there are strong devia- 
tions between theory and experiment. 

The simplest explanation for the smaller weight in the 
experiment is that some of the structures that contribute 
at high energy may be broad and assigned to the back- 
ground (see Ref. l3l|) . We should say that because we 
neglect life-time effects even the sharp structures of our 
computation will be much broader in reality and this will 
be more important as the energy increases since, quite 
generally, the phase space for decay processes increases. 
Also high-energy excitations will be much more sensitive 
to disorder since they involve short wave-lengths and this 
will also tend to make them broader. The broad struc- 
tures will be assigned to the background explaining the 
difference between theory and experiment. We should 
also take into account that can be more depressed 
due to the larger impact of quantum fluctuations in the 
metallic phase. 

In order to have a broader view of the trends in the evo- 
lution of the spectral weight, we show in Fig. 1131 a com- 
pilation of measurements of S'^^-^ from different groups 
and at different dopings. For x — 0,0.125 and 0.14 the 
overall weight seem to be roughly conserved whereas for 
X = 0.16 the detected weight seems to be significantly 
smaller. It will be interesting to see weather this is due 
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to a transfer to high-energies where absolute measure- 
ments are not available. 



V. CONCLUSIONS 

In this work we have analyzed the distribution of spec- 
tral weight in neutron scattering experiments in cuprates. 

In the insulator the spectral weight has been estimated 
in detail in the Heisenberg model and then translated into 
a Hubbard response with the use of an effective shield- 
ing factor. We find an effective shielding factor which 
essentially coincides with a proposal of Hirsch and Tang 
obtained by analyzing numerical data^ In addition we 
have estimated, with the aid of the TDGA, the weight at 
the energy of the charge transfer transitions. 

The inelastic longitudinal spectral weight at low ener- 
gies is in broad features and thus hard to detect, although 
progress has recently been made^i The transverse part 
has the well know quantum renormalized propagating 
spin wave modes which will be the features more accessi- 
ble to experiment. Theory and experiment agree well in 
the insulator. 

For the doped phase we computed the dynamical re- 
sponse within the TDGA. We considered metallic stripes 
but we expect our results for overall weights to be largely 
insensitive to the specific texture. One finds extra broad 
features due to the particle-hole continua plus propagat- 
ing collective modes. In analogy with the insulator we ex- 
pect the latter to be reduced in intensity due to quantum 
fluctuation and to dominate the spectral weight reported 
in experiments. 



Our motivation was to understand why such a small 
fraction of the naively expected spectral weight is actu- 
ally detected. We see that after taking into account all 
the reduction factors we passed from a weight which was 
an order of magnitude smaller than expected to a weight 
which is ^ 2/3 of that expected. Given the uncertainties 
involved we think this value is reasonable. 

By far the most important factor in reducing the de- 
tectable weight turns out to be the fact that a large frac- 
tion of the spectral weight is in broad features or at too 
high energies and this make the detection of the total 
spectral weight experimentally very challenging. In the 
doped phase a consistent fraction of the transverse spec- 
tral weight is predicted to lay at intermediate energies 
which may become accessible to experiment in a near 
future. 

We see no reason why the broad/high energy modes 
can not have a strong impact on the effective attraction 
between holes relevant for superconductivity and there- 
fore their experimental characterization remains an im- 
portant open problem. 
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